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ABSTRACT 

We study the internal color properties of a morphologically selected sample of spheroidal galaxies taken from 
the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) ERO program of UGC 10214 ("The 
Tadpole"). By taking advantage of the unprecedented high resolution of the ACS in this very deep dataset we 
are able to characterize spheroids at sub-arcseconds scales. Using the Vgo6 and 7si4 bands, we construct V — I 
color maps and extract color gradients for a sample of spheroids at I&\4w < 24 mag. We assess the ability of 
ACS to make resolved color studies of galaxies by comparing it with the multicolor data from the Hubble Deep 
Fields (HDFs). Here, we report that with ACS/WFC data of £ 10x less exposure time than in the WFPC2 
HDFs it is possible to confidently carry out resolved studies of faint galaxies at similar magnitude limits. We 
also investigate the existence of a population of morphologically classified spheroids which show extreme 
variation in their internal color properties similar to the ones reported in the HDFs. These are displayed as 
blue cores and inverse color gradients with respect to those accounted from metallicity variations. Following 
the same analysis we find a similar fraction of early-type systems (~ 30% — 40%) that show non-homologous 
internal colors, suggestive of recent star formation activity. We present two statistics to quantify the internal 
color variation in galaxies and for tracing blue cores, from which we estimate the fraction of non-homogeneous 
to homogeneous internal colors as a function of redshift up to z < 1 .2. We find that it can be described as about 
constant as a function of redshift, with a small increase with redshift for the fraction of spheroids that present 
strong color dispersions. The implications of a constant fraction at all redshifts suggests the existence of a 
relatively permanent population of evolving spheroids up to z < 1 . We discuss the implications of this in the 
context of spheroidal formation. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: fundamental parameters — galaxies: 
evolution 



1. introduction 

In recent years, the use of field elliptical galaxies has be- 
come a conventional tool for testing between very different 
structures of formation, where the main competing views 
of galaxy formation at high and low redshift are often re- 
ferred as monolithic or early formation and hierarchical or 
late formation (Peebles 2002). Many observational stud- 
ies have been devoted to study a number of "scaling rela- 
tions" in ellipticals i n rich clusters; i.e. the low scatter in the 
fundamental plane jb iorgovski & Davis|l2£7J tPrgs sler et al.l 
119871 iTreu et alJ Il999t IvanD okkum & Franxlll996l) and in 
the color-magnitude relation (CMR) (Bower et al.l 119921 
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Ellis etai]fl997t landage & Visvanathar3 ll978l) . These im- 
ply a high degree of homogeneity in the stellar population 
which reinforced the idea of a monolithic collapse model in 
which elliptical s formed during a rapid collapse at high red- 
shift JE ggen et al.ll962l) . 

The simple view of early formation contrasts with the pre- 
dictions of models where galaxies assemble hierarchically as 
the result of the merging of smaller sub-units, a t a rate gov- 
erned by the merger of cold dark matter halos (Baugh et alJ 
119961 IWhite & Frenklll991t IWhite & Reeslll978l) . Although 
the prevailing view in the past was that ellipticals in the field 
formed in isolation in a high-redshift initial burst of forma- 
tion, as did their clustered counterparts, several authors have 
shown that the observational properties of field spheroids 
(collectively E and SOs) are incompatib le with the bulk of 
the population forming at high red shift (Barg er et al.lll999t 
Menantea ul-tafll 1 999t iTrager et al Jl2000t IZepfll 1 997lh Most 
of these studies were prompted by the success of hierarchi- 
cal models in the predictions of a large quantity of observ- 
able properties of galaxies, suggesting that massive field el- 
lipticals could only have been assembled recently (z < 1). 
The general consensus that arose is that a single short pe- 
riod of formation at high redshift is incompatible with the 
observations. There is growing evidence of the existence of 
blue spheroids with perturbed colors at intermediate redshift 
and a deficit of red systems compared to monohthic mod- 
els predictions (iK odama et al.|| 199 91 iMenanteau et al JI200 1 al 
ITreu & Stiavellilll999t IZepfll 19971) More recently, however, 
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IB ell et ai] ( 120041) have reported using a large sample of HST 
galaxies that E/SOs dominate the red-sequence of galaxies at 
Z ~ 0.7. As the outcome of these new studies a new view of 
elliptical formation is emerging. 

A surprising result from this wave of studies is the exis- 
tence of field spheroids with blue central regions at interme- 
diate redshifts 0.4 < z < 1.0. These objects were orig inally 
detec t ed from their V^pf, — 7rh color maps (Abraham et al. 
119991 iMenanteau. Abraham. & Ellisl 1200 lU) in the Hubble 
Deep Fields (HDFs), and are suggestive of recent star- 
formation activity associated to central region of the galaxy. 
The presence of blue core spheroids in the HDFs has gen- 
erated interest in reproducing t heir observed properties us- 
ing both semi-analytic models ( Benson, Ellis, & Menanteau 
1 20021) and extended monolithic collapse ([Jimenez et al J19981 
Menanteau, Jimenez, & Matteucci 200 f3). The hie rarchi- 
cal description of Benson, Ellis, & Menanteau (2002) accu- 
rately predicts the number of spheroids and the degree of 
color variations observed in spheroids, but misses most of 
the red spheroids in the sample. On the other hand multi- 
zone chemo-dynamical models can predict the existence of 
blue cores and homogeneous colors of evolved ellipticals 
with success, by adjusting the redshift of formation of the 
galaxy, although the cosmological mechanisms responsible 
for detailing the f ormation of spher oids are somewhat un clear. 
jFri^c a & Terlevich1l200 It [Menanteau. Jimenez. & Matteuccil 
12001 hi) 

The use of resolved colors in evolutionary studies, although 
largel y unexploited f or characterizi ng dis t ant galaxies (see 
Abrah am etallll999t iHinklev & Iml 120011 IMenanteau et all 
2001a; Tamura et al. 2000| for examples) supplies a new set 
of independent constraints to the traditional use of counts and 
redshift distributions alone, particularly useful when attempt- 
ing to discriminate between models of elliptical formation. 
Moreover, it is appealing to examine independent spheroid 
samples that look for galaxies with color inhomogeneities 
like those reported in the HDFs. Over larger datasets, they 
can provide important clues for determining the fraction of 
spheroids still experiencing star formation. Unfortunately, un- 
til the arrival of the ACS, resolved colors analysis have been 
limited only to superb signal-to-noise and highly time con- 
suming datasets such as the HDFs. In this paper we explore 
the advantages arising from exploiting the unprecedented ca- 
pabilities of the ACS in obtaining HDF-like datasets with only 
a fraction of the integration time. We use the first data avail- 
able from the ACS to study the resolved color properties of 
faint distant galaxies and compare them with previous results 
from the HDFs. 

A plan for the paper follows. In Section 2, we describe 
the observations and data reduction of the ACS images and 
our morphological classification of galaxies utilized in segre- 
gating spheroids. In Section 3 we discuss the resolved color 
properties and the construction of color maps of spheroids in 
our sample. In Section 4 we present our methodology for 
characterizing the resolved color properties using a model- 
independent approach. In Section 5, we discuss the results of 
our analysis and finally in section 6 we summarize our con- 
clusions. 

2. SPHEROIDS IN THE ACS/WFC UGC 10214 FIELD 

The spheroids in our study were selected from one of the 
first s cience observations acquired by the ACS (tFord et alJ 
12001 Wide Field Channel (WFC) as part of the Early Re- 
lease Observations (EROs) program. The ACS/WFC houses 



two 2048 x 4096 butted CCDs separated by a ~ 2" gap pro- 
viding a FOV of 202" x 202" with a pixel size of 0.05". 

2.1. ACS/WFC Observations 

In our analysis we utilize the deep images of UGC10214 
(also know as "The Tadpole", VV029 and Arp 188), a bright 
spiral with long tidal tail at z = 0.03136, imaged in the 
F475W(g), F606W(V) and F814W(7) filters. Due to an er- 
ror in telescope pointing observation in which the "head" of 
the Tadpole was cut off the ACS field of view, a second set 
of observations was performed soon after, leading to an ex- 
tra deep exposure on the overlapping region. In Tabled we 
report the observations of UGC 10214, including exposure 
times, number of orbits and depth. The final observations re- 
sulted in a superb multicolor dataset in which the depth limit 
is <1 mag shallower than the southern HDF (HDF-S) for all 
the observed bands. The properties of the Tadpole galaxy it- 
self and its formation act ivity in young st ars in the tail have 
been studied in detail by Tran et al l (120031) . In this paper, we 
focus instead on the properties of a sample of galaxies in the 
Tadpole field. It is important to point out that although the ap- 
parent large size of the Tadpole galaxy fills the whole frame, 
there is a large remaining area of the WFC image of 11.54 
arcmin 2 which is essentially clean of any foreground interfer- 
ence and hence perfectly suited for the study of distant field 
galaxies in a similar fashion to those performed in the HDFs. 
Hereafter we will concentrate on and refer only to the objects 
in the background of the WFC image of UGC 10214. In fig- 
ure [0 we show a composite color image of the field in which 
it is possible to appreciate the abundance of faint background 
galaxies. 

2.2. Image Reduction and Object detection 

The d ata were fir st processed by the STScI CALACS 
pipeline (Hack 1999), which included bias and dark subtrac- 
tion, flat fielding and counts-to-electrons conversion. The 
images were then processed by Apsis (ACS pipeline sci- 
ence investigation software) at Johns Hopkins University. 
This pipeline which measured rotations and offsets between 
dithered images from both pointings, rejected cosmic rays 
(CRs) and combined them into single geometrically corrected 
images using the drizzle method in each band and a combined 
image used for obje ct detection. The reader is referred to 
Blakeslee et all fc003l) for a thorough description of Apsis. 

Object detection, extraction and integrated p hotometry 
were taken from the SExtractor (Bertin & Arnouts 1996) cat- 
alogs produced by the Apsis pipeline. Extensive simulations 
and tests were conducted to determine the optimal parame- 
ters for extraction and deblending of ACS/WFC field galaxy 
sources. A detailed description of th e procedure and parame- 
ters is given by Bemtez et al. (2004). In a nutshell, object ex- 
traction was carried out over a specially created detection im- 
age composed of the addition of the images in all bands. Ob- 
jects detected above a certain threshold limit within the detec- 
tion image are included and selected for subsequent photom- 
etry. Apsis provides -through SExtractor- a stream of magni- 
tude measurements computed for each of the bandpasses. For 
the purpose of our analysis and object selection we choose the 
near-total magnitude MAG_AUTO based on Rron-like ellipti- 
cal aperture magnitude measured within k times the isophotal 
radius. While the SExtractor magnitudes are calibrated in the 
AB system as Apsis products, we transformed them into Vega 
magnitudes to ease the comparison with previous works, em- 
ploying the photometric zero-points from the ACS photomet- 
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Table 1. Description of UGC 10214 observations 













N orbits Area (arcmin 2 ) " 




Image depth 




Filter 


RA(J2000) 


DEC(J2000) 


Expt time(s) 


N exp 


Overlapping 


Shallow 


HDF-S 


F475W 


16:06:17.4 


55:26:46 


13600 


12 


6 11.54(14.48) 


27.64 (27.86) 


27.27 (27.49) 


27.97 (27.90) 


F606W 


16:06:17.4 


55:26:46 


8040 


12 


4 11.54(14.49) 


27.48 (27.70) 


27.14(27.36) 


28.47 (28.40) 


F814W 


16:06:17.4 


55:26:46 


8180 


12 


4 11.54(14.46) 


27.03 (27.25) 


26.62 (26.84) 


27.84 (27.77) 



Note. — "The values given in this column are the effective (i.e., "Tadpole-less") areas used in the analysis, and the accompanying values in parentheses are the total areas actually 
covered by the ACS observations. fc HST image depths taken from Benftez et al. ' 2004) for point-like and extended (values in parentheses) sources. 



Fig 1 available at http://acs.pha.jhu.edu/~felipe/e-prints/Tadpole 



FIG. 1. — Composite color image of UGC 10214 ("The Tadpole"). The image illustrates the abundance of background objects observed in this deep ACS 
image. The impressive depth and image resolution is shown in the figure insert. 



ric calibration program ( Siri anni et alJl2004h . Hereafter we 
will refer to magnitudes in the Vega system. 

2.3. Selecting Spheroids 

In selecting field spheroids we first choose objects with 
integrated magnitude I%\4 < 24, using MAG_AUTO from 
the SExtractor products. This is our starting point for 



segregating early-types. We decided on selecting objects 
bright er than I%u < 24 — the sam e criteria adopted by 
Menanteau, Abraham, & Ellis (2001a) in the HDFs — as this 
enables both reliable morphological classification and abun- 
dant signal-to-noise pixel information on deep HST images 
and will provide us with the right tools to compare with the 
HDFs sample later. This lead us to an initial sample of 373 
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Table 2. ACS Photometric Zero Points in the AB and Vega systems 



Filter 


AB System 


Vega 


F475W 


26.043 


26.144 


F606W 


26.505 


26.410 


F814W 


25.941 


25.496 



Note. — ACS/WFC photometric zero-points obtained from Sirianni et al. ' 2004). 



objects. 

We employ th e same procedure for selecting sphe roids 
as described in | Menanteau7A braham. & Ellis (2001a|) and 
Abra ham et alJ (1 19961) based on a combination of visual clas- 
sification and machine morphological analysis. Stars were 
initially removed using the SExtractor star class parameter 
in addition to visual inspection of all objects in the sample. 
Here, we briefly summarize the strategy taken and refer the 
reader to these works for a full description of the methodol- 
ogy used. Galaxies in the Tadpole field were morphologically 
classified using an automated classification based on both the 
central concentration ( Cj and a symmetry (A ) parameters (A - 
Q from Abrah arrietall Jl996), as well as using visual classi- 
fications made by one of us (FM). Visual classification have 
been shown to agre e very well with A - C classes in previ - 
ous studies ( Brinchmann et al. 1998; Menanteau et al. 1999). 
The segregation of early-type systems is particularly robust as 
their chief diagnostic parameter for discrimination is central 
concentration. It is worth noting the key advantage in adding 
A - C to visual classes alone, as they represent objective mea- 
surements of the morphological properties of galaxies and can 
be modeled and easily reproduced in the future if required. In 
building our final sample of spheroids we used a combined 
criteria of using both A - C and visual classes, which we be- 
lieve is a robust approach. In figure we compare the A-C 
values computed for galaxies keyed to their visual classes in 
three broad categories: E/S0, spiral and irregular. We selected 
spheroid galaxies as objects which were visually classified as 
E, E/S0 and SO in the Hubble scheme. When there were dis- 
crepancies between visual classes and A - C, they were exam- 
ined in detail, to ensure an unbiased selection. In this fashion 
we constructed a final catalog containing 116 ACS spheroids, 
which is ^ 43% larger than the only other HST existing cata- 
log; the HDFs fields with 79 spheroids. 

2.4. Bayesian photometric redshifts 

Despite the absence of spectroscopic information for galax- 
ies in the Tadpole field, we enhance the sample with the pho- 
tometric redshift estimates computed using the Bayesian Pho- 
tometric Redshift package (BPZ) (Bemtez 2000) included in 
the Apsis pipeline products. BPZ follows a Bayesian statisti- 
cal approach to estimate redshifts employing a set of galaxy 
templates and magnitude-redshift priors. While it is often 
assumed that several filters are necessary to achieve accu- 
rate measurements, the addition of priors can lead to reli- 
able redshifts with only the F475W, F60 6W and F814W fil- 
ters. As described in Bemtez et al. (2004), when the Bayesian 
methodology is calibrated with the Northern HDF (HDF-N) 
spectroscopic sample, it has an excellent behavior at z < 1.5, 
without catastrophic outliers and much better performance 
th an maximum likel ihood alone. Furthermore, according 
to Bemtez et al. (2004), simulations of the efficiency of the 
photometric redshifts as a function of magnitude confirmed 
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-0.8 -0.6 -0.4 -0.2 

log C 

FIG. 2. — The computed values of the Central Concentration and Asymme- 
try and their visual classes. Plot symbols denote E/S0 systems as red ellipses, 
spirals as green spirals and irregulars as blue asterisks. 



that for bright objects, I$\4 <24 the tipical uncertainties are 
dz/(l +z) ~ 0.1, which coincides with the magnitude limit 
set for our sample where photometric estimates are most re- 
liable. Despite the small estimated errors, we opt to take the 
safer approach of using the photometric redshift only over the 
integrated properties of the sample rather than those of indi- 
vidual galaxies. 

3. RESOLVING SPHEROIDS 

Understanding the behavior of the point-spread function 
(PSF) becomes important when studying distant objects at 
pixel-scales. In the present analysis our main concern is 
that PSF variations as a function of wavelength may lead to 
spurious centrally concentrated inhomogeneities arising from 
some centrally concentrated profiles in some spheroids. To 
study this effect we made extensive simulations of the effect 
of the PSF by artificially creating spheroids which were sub- 
sequently analyzed using the same method applied to the ob- 
served data. These are presented in AppendixlAl We conclude 
that the effect of the PSF will not significantly influence any 
of the indicators used in our analysis. 

Prompted by th e findings of spheroids with ce ntral blue 
cores in the HDFs (Menanteau. Abraham. & Ellis 2001a), the 
first step in our analysis is to construct Vgo6 —^814 color maps 
for all the ellipticals in our sample. This color peculiarity has 
only been reported for early-types in the HDFs so far, and it 
is compelling to verify whether it will also be present in other 
unrelated samples such us the current one. For this, we com- 
pute pixel maps, surface brightness, jj., for all galaxies in all 
three bands, using the well-known transformation, 

jll = zp — 2.51og(counts/exptime) +51og(pixelscale) 

where zp is the zero-point for a given filter in the Vega system 
from Table|2 the pixelscale is 0.05 arc/pixel and counts is the 
number of electrons in each pixel. Subsequently, color maps 
were computed, selecting only contiguous pixels above a sur- 
face brightness threshold of /I606 < 25 mag/arc sec 2 . We note 
that we use this surface b right ness limit only in the computa- 
tion of 8 (V — /) in section RTTl in order to maintain consistency 
whe n co mparing with previous observations. However in sec- 
tion we w iH discusss a more robust method for selecting 
and measuring the light of spheroids. 

3.1. Old ellipticals vs Blue ellipticals 
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FIG. 3. — Compaiison of the V606 ~ ftl4 internal color properties of 
spheroids in the Tadpole field and in the HDF-N. The top two galaxies are 
taken from the HDF-N sample and the ones below, from the ACS sample. 
The left column shows examples of blue nucleated spheroids and on the right 
normal spheroids. 



Our first objective is finding a population of spheroids with 
blue cores and inhomogeneous internal colors in the current 
ACS sample mirroring those found in the HDFs. After an 
initial inspection, we conclude that in most cases spheroids 
show fairly smooth and homogeneous Veo6 — hu color maps, 
but the existence of a subsample with perturbed internal col- 
ors is quite evident. We report a fraction of spheroids with 
blue cores comparable to the HDFs in our ACS sample, with 
similar sizes and strength in their color differences. It is im- 
portant to remark that this is the first sample, apart from the 
HDFs in which resolved studies of distant galaxies are possi- 
ble. To illustrate the similarities in our findings, in Fig. [3] we 
show examples of spheroids with normal and peculiar internal 
color properties in the HDF-N and the present sample of ACS 
spheroids. The similarities in the form of the blue nucleated 
ones in both samples are quite evident. 

4. METHODOLOGY AND ANALYSIS 

After considering the color maps of spheroids, we want to 
investigate how the internal colors of field ellipticals vary as 
a function of redshift and quantify the fraction that departs 
from a passive evolving population. Although we can iden- 
tify a distinctive population of spheroids with blue cores and 
color inhomogeneities from visual inspections alone, we wish 
to develop a method for quantifying its strength in an ob- 
jective way. A standard tool for this purpose is to model 
thei r color p roperties using stellar synthesis population mod- 
els (|Abrahamet]d][l999l [MinintiMlIiniiOOlB- Unfortu- 
nately, these predictions are quite sensitive to the galaxy red- 



shift of observation and our present sample only encompass 
photometric estimates, making this approach toward charac- 
terizing individual galaxies somewhat uncertain. Instead we 
chose to explore a model-free approach and focus on the 
integrated properties, which in our view represent a safer 
way to evaluate the trend of the population. For this pur- 
pose we revisit the internal color dis persion estimator from 
Menanteau, Abraham, & Ellis (2001a) and devise a new mea- 
surement based on the slope of the color gradients of a galaxy 
that we will show is more sensitive to the presence of central 
blue spots. 

4. 1 . The internal color dispersion 8(V — I) 

The idea of studying the internal colors of an individual 
galaxy can be considered as a generalization of the use of 
the color-magnitude relation i n the study of the history o f 
elliptical galaxies in clusters ( Bower. Lucev. & Ellis 1992), 
where the photometric dispersion is employed as a power- 
ful mechanism to determine variations in their star formation 
history. The first quantity that we present is the signal-to- 
noise weighted scatter in the inter nal V^nts — hu colors of a 
galaxy, 8(V — I), as defined in lMenanteau. Abraham. & Ellis! 
(2001a). This measurement has proven a good tracer of pe- 
culiarities in spheroid colors at pixel scales in the sense that 
objects with low values of 8 (V — /) have very uniform and 
normal color maps, while those with higher values do present 
disturbed internal colors and in many cases blue cores. 

A key motivation for re-visiting this estimator is that it has 
already been applied and calibrated with a similarly selected 
sample of galaxies in the HDFs, where it proved to be a good 
discriminant of objects with blue cores vs ellipticals with nor- 
mal colors. This feature enable us to make direct comparison 
between the HDFs sample 8(V — I) measurements and the 
current one. In addition, it is a model independent method for 
studying differences in the star-formation histories of galax- 
ies, making it free of assumptions regarding the epoch of for- 
mation in ellipticals. 

In our estimates of 8{V — I), we follow the same proce- 
dure used in the HDFs, and we also chose to concentrate 
only on the V606 and Ig 14 bands for computing color dis- 
persions, as the resulting V(,Q6 — hu color contains signif- 
icantly smaller observational errors than #450 — V606 , spe- 
cially for systems dominated by old stellar populations (see 
Menanteau. Abraham. & Ellis 2001a, for details). However, 
it is important to note that strong variations in the internal 
colors are also present when looking at the the #450 — Vso6 
colors. We select all contiguous pixels within a surface bright- 
ness limit of < 25 mag/arcsec 2 in V6O6, to isolate the galaxy 
from the background sky and maximize the S/N ratio per pixel 
associated. Using only the pixels within this limit, 8 (V — I) 
is computed from the weighted signal-to-noise distribution of 
Vso6 — hu colors for individual pixels, which in turn identi- 
fies the internal homogeneity of a galaxy, using the following 
recipe: 

8(V-I) ^ Xi ~f±l SNRixi) (1) 



S>,-S(s,-)SNR(s f ) 
ES(x,)SNR(jc,) 



(2) 



where S(x) is a selection function for pixels with signal/noise 
above a certain threshold such that S(x) = 1 for pixels above 
and S(x) = for pixel below the threshold. SNR(xi) is the 
signal-to-noise ratio for a given pixel color Xj. The selection 
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function S(x) and the weighting according to SNR(x) address 
biases arising from noise variations at pixel-scales by reject- 
ing low signal pixels and weighting the pixels' contribution 
proportionally to their signal. 

4.2. The mean {V(,Q6 — hli) ( r ) slope of the galaxy 

In addition to the color dispersion we wish to develop an al- 
ternative quantity to isolate spheroids with non-homogeneous 
internal colors. This should be suited to measure the strength 
of the blue cores and not prone to dispersions arising from par- 
ticularly reddened core. For this we concentrate on a method 
that quantifies the variation in the color as a function of the 
galaxy radius. Given that in most of the cases, color pecu- 
liarities manifest themselves via the presence of blue cores, a 
clear pattern of this effect is the display of inverted (i.e. pos- 
itive sign) color gradient, as opposed to "normal" passively 
evolving ellipticals which either show flat or slightly negative 
color gradients, as expected from metal enriched cores. 

We focus on the V(,Q(, — hi 4 colors to compute color gra- 
dients as a function of the radius r, as these represent the 
bands where the signal-to-noise is higher and can be more 
directly related to 8(V — I). In the extraction of the color gra- 
dients, initially we measure the centroid and second order mo- 
ments (x 2 , y 2 , xy) of the galaxy utilizing the 7si4 band, from 
which we determine the ellipticity parameters of the objects 
(a,b,9). It is worth noting that we keep the same elonga- 
tion ratio e = a/b up to the isophotal limit to which the color 
gradient is computed. Next, using concentric ellipses, we cal- 
culate the V<506 — hu( r i) gradients as the median color in the 
shell between r,_i and r i+ i. To gauge the slope of the color 
gradient, we fit a fourth order polynomial function to the ob- 
served (V606 — hu)( r )> U P to a certain maximum radius r maj: , 
from where we measure the slope by simply obtaining the first 
derivative of the fitted polynomial function. In order to con- 
sistenly measure the slope of galaxies over similar physical 
areas on spheroids and at different redshifts, we compute the 
slope up to a rad ius r mar following a si mila r aproach to the 
one de scribed by Papovi ch et al.l (|2003) and Conseli ce et al.l 
(2000). This consists of choosing the aperture radius r max of 
the galaxy as r,„ nx = 1.5 x r p with (tj = 0.2) where r p is the 
iPetrosianH 19761) radius and rj(r) is defined as: 



rj(r)=/(r)/(/(r)> 



(3) 



where I(r) is the surface brightness of the galaxy in an annu- 
lus of radius r and (I(r)) is its mean value up to the same ra- 
dius. The Petrosian radius depends only on the galaxy surface 
brightness and it is independent of the redshift of observation. 
Formally, the mean (V606 — ^814) (>*) slope of the galaxy can be 
expressed as: 



slope = 



Li-(V-I)(r)\n-&n 



(4) 



where the derivative is taken over the fitted function, Ar, is 
the size of the concentric shells and r is measured in kpc. We 
repeat this procedure to all the galaxies in our sample. In 
Section|5]we will investigate how it performs compared with 
8 (V — I) and present examples. 

4.3. Comparing the indicators 

An important objective of our analysis, is to carry out di- 
rect comparisons of the internal color dispersions between the 
present ACS sample and those observed in the HDFs. For this 




^(^606 I B14 ) 

FIG. 4. — Comparison between the two statistics used to probe varia- 
tions in the internal color properties, S(V~ I) and the mean slope (V — /) 
for spheroids in the Tadpole field, represented by open circles. Values for the 
HDF-N are also shown as hatched circles. The size of the symbols is keyed 
to the Ig 14 magnitude as indicated in the panel within. Galaxies with large 
positive mean slope (V^ — h\4)(r) bear blue cores, while those with large 
S (V — I) values show strong color dispersion in general associated with blue 



purpose, we will include in our study the morphologically se- 
lected catalog of spheroids in HDF-N from Menantea u et alJ 
(2001a). The HDF-N sample was constructed using the same 
selection limits and morphological classification, ensuring a 
robust comparison between samples. We decide to include 
only the Northern HDF as this contains the largest num- 
ber of spectroscopic redshifts, making possible an investiga- 
tion of how the internal color dispersions vary as a function 
of redshift as opposed to the photometric estimates in ACS 
spheroids. After computing 8(V — I) and the mean slope for 
all galaxies in the ACS and HDF-N samples we examine how 
they perform against each other and subsequently compare 
with HDF-N spheroids. 

Figure |4] shows the values of 8 (V — I) versus the mean 
slope (V606 — h\4)(r) for all spheroids in the Tadpole field. 
We can observe a clear correlation between both statistics for 
the larger part of the sample: galaxies with high positive mean 
slopes (indicative of blue cores) also show high values of 
8(V — I) (i.e. high internal color dispersions). Although the 
inverse relation is true for most galaxies (i.e high color disper- 
sion values (8(V — /)>0.06) also have high positive slopes) 
we see in Fig.|4]that several high 8(V — I) galaxies have flat 
or even negative slopes in some cases. This results as galaxies 
with metal enriched cores tend to have redder centers, mak- 
ing their mean slopes slightly negative or flat, while keeping 
higher than normal values of 8(V — I). From Fig.|4l we wish 
to highlight that the mean slope (V606 — /si4)(r) is a more ro- 
bust statistic for tracing blue cores as spheroids with positive 
slopes and hence blue cores exhibit high internal color dis- 
persion. Also in Fig. we present the values computed for 
the HDF-N spheroids. It is reassuring that ellipticals populate 
the same region in the 8 (V — I) — median slope space for both 
samples. After individually inspecting the objects with high 
8 (V — I) values but low or negative slopes, we confirm that 
these correspond to objects with reddened cores. We conclude 
that although both statistics succeed in tracing the observed 
color variations, the mean slope is a more robust statistic for 
tracing blue light in spheroids and is not prone to color differ- 
ences arising from color gradients normally accounted for in 
passively evolved ellipticals. 
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Redshift 

FIG. 5. — The scatter S(V — I) as a function of redshift for spheroids in 
the Tadpole field are shown as open circles. For reference the values for 
spectroscopic redshifts from the HDF-N are also shown as filled triangles. 
Horizontal bars at the top represent typical 1 o~ photometric redshift errors for 
galaxies in the ACS sample. 



5. RESULTS 

Establishing the existence of a population of field spheroids 
with marked color dispersions, resembling those exhibited by 
HDFs spheroids is a major objective of this paper. From the 
constructed color maps alone we can conclude that this is a 
property shared by field spheroids at intermediate redshifts 
regardless of the sample in study. The presence of strong 
variations in the internal color distributions suggests episodes 
of recent start-formation activity. However, detailed measure- 
ments regarding the intensity and timescales involved requires 
precise spectroscopic information, not available for the sam- 
ple in our study. Instead, we opt to determine in a model- 
independent way the proportion of active systems as a func- 
tion of redshift and based on this information, make educated 
guesses regarding the evolution of early-type galaxies. 

Initially, we wish to investigate the behavior of S(V — I) 
as a function of redshift. For this, in figure [5] we display 
the computed values of 8 (V — I) as a function of their pho- 
tometric estimations for the Tadpole field (open circles) and 
spectroscopic redshifts for the HDF-N (filled triangles). Error 
bars represent Iff estimates obtain ed using the bootstrap re - 
sampling technique introduced by Efron & Tibshirani ( 1986). 
From figure[5]we observe how the fraction of spheroids with 
high 8 (V — I) values slightly increases as a function of red- 
shift in the same way as reported for HDF-N spheroids. We 
make a quantitative measurement and focus on the fraction of 
systems with high 8 (V — 7)values compared to low ones as 
this measurement is not prone to selection effects due to the 
redshift distribution of spheroids in our sample. We separate 
the sample into two redshift bins, such as < ZbM < 0.6 and 
0.6 < Zbini < 1-2 and calculate the fraction of spheroids with 
8(V — I) > 0.09 to the total number. We report an increase 
from 28% in ZbM to 39% in Zbirii- However, we point out 
that this might also be consistent with no increase given the 
relatively small number of objects involved. 

It is reassuring to confirm that both samples show a very 
similar trend in the observed spread as well as its increase with 
redshift. Assuming that in most cases high values of 8(V — I) 
are associated with recent star formation activity within the 
galaxy, the rise of 8 (V — I) with redshift indicates a clear 
trend in this sense. For the HDFs, [Menanteau et alJ {2001a) 
have computed simple models that can explain the existence 



of galaxies with high 8 (V — I) as the result of small additional 
bursts of star formation activity. 

If we assume that strong internal color dispersion are indeed 
associated with recent star formation, it is interesting to ob- 
tain the fraction of evolving systems as a function of redshift. 
As observed from figure |4] the mean slope (V606—hu)( r ) 
is better suited to probe blue cores and internal inhomo- 
geneities than 8(V — I), therefore we choose to use this in- 
stead of 8 (V — I) for estimating the ratio of spheroids with 
strong internal color variations to the total number of ob- 
jects. In figure[6]we show computed values of the mean slope 
(^606 — hu){r) as a function of redshift for both the Tadpole 
field and the HDF-N sample, using the same symbols as in 
figure[5] For completeness we also display the values for the 
HDF-N, as this allows to make consistency checks for vari- 
ations that could arise between samples, in part due to the 
use of photometric redshifts. Error bars show Iff computa- 
tions calculated by carrying out Monte Carlo simulations over 
the data. To illustrate how the mean slope is used in probing 
galaxies with blue cores, in figure[6]we also show two galax- 
ies of contrasting internal colors and the very distinctive val- 
ues of the mean slope they have. It is interesting to note from 
this figure how galaxies tend to populate rather uniformly the 
redshift space. 

5.1. The mean slope (V<506 — ^8l4)( r ) redshift dependence 

In order to reliably make use of the mean slope 
(^606 — hu){ r ) m tracing blue light in spheroids, we need to 
investigate its variation as a function of redshift introduced 
from measuring the mean slope using the Va)6 and /sh filters. 
We attempt to reproduce through simulations the observed 
spatially resolved color properties of spheroids with blue nu- 
clei. Our key assumption when modeling these spheroids is 
that at intermediate redshifts they can be described using a 
two component model: a young stellar central component 
(responsible for the blue light) superimposed on an old sys- 
tem. We modeled the integrated V(,06 — hu colors of both 
stellar component using spectral energy distributions (SED) 
from the Bruzual & Chariot (1996) (BC96) library. For the 
old stellar component we chose an SED corresponding to ex- 
ponentially declining star formation (e-folding time 1 Gyr) at 
age 12 Gyr, which resembles the colors of a normal ellipti- 
cal at z ~ 0. For the central region we employ the SED of 
a recent star-forming elliptical as observed after 1 Gyr since 
its formation. For each component their observed V606 — /gi4 
colors and apparent magnitudes are computed as as a func- 
tion of redshift (^-correction only) assuming a flat universe 
(H = 71 km s-'Mpc" 1 , Q. Q = 0.27, Q. A = 0.73). 

To emulate the geometrical properties of the modeled 
galaxies (as those of the blue nucleated spheroid in Fig. |6j 
we created artificial ellipticals from the two components using 
a customized version of the IRAF package ARTDATA. We 
used a fixed physical size of half-light radius r e = 3 kpc for 
a de Vaucouleurs profile for the old component. We assume 
that the blue central component corresponds to only 15% of 
the total stellar mass of the galaxy M — 3 x 10 M©, as this 
supplies enough blue light to reproduce the observed blue col- 
ors, and effective radius of 0.15r e . 

For the resulting artificial galaxies we compute the mean 
slope (V606 — hi4){r) between 0.2 < z < 1.2 using the same 
procedure as for the real data. We show the results of the exer- 
cise in figure[6] The solid line represents the computed values 
of the mean slope (Ve06 — ^814) M recovered from the simula- 
tions. We observe that for simulated spheroids the mean slope 



8 



F. Menanteau & The ACS Science Team 






FIG. 6. — The mean slope (Vgos - hu)^) as a function of redshift for spheroids in the Tadpole field are shown as open circles. For reference the values 
with spectroscopic redshifts from the HDF-N are also shown as filled triangles. The dotted line is the limit used for selecting evolving systems as a function 
of redshift. To illustrate how the mean slope (V606 — ^814)('°) traces core colors, we show two examples of objects with positive/high (top) and negative slopes 
(bottom), with their associated position in the mean slope vs redshift space. We also display their Vf,o6 ~~ '814 color maps and their corresponding (Vgo6 — hu)( r ) 
extracted gradients in the far right panels, with their associated Iff error bars. The solid red line in the extracted color gradient represents the fitted slope used in 
the computation of the mean slope. The solid line shows the mean slope values obtained for simulated galaxies as a function of redshift. 



changes little as a function of redshift, with a bump at z ~ 0.65 
which coincides with the peak in the A(V606 — h\4) colors dif- 
ference between both components. We conclude that although 
our filter set is more favorable to trace blue light near z ~ 0.65, 
the variations as a function of redshift are small enough that 
this will not significantly affect our results. 

5.2. The fraction of active systems 

When computing the active to total galaxies ratio, we need 
to calibrate the mean slope (V606 — hu)(f) to differentiate 
between active versus passively evolved systems. Certainly 
high values of the mean slope are directly linked to large in- 
ternal color variations, and for the high mean slope vs. red- 
shift envelope, is where undoubtedly the most striking cases 
of blue cores arise. On the other hand, galaxies with small 
values, near zero or negative show no significant color vari- 
ations. However, defining a limit that separate active vs pas- 
sive systems is a rather uncertain and imprecise exercise. To 
assess this problem, we opt for interactively calibrating this 
limit, by visually inspecting all galaxies and their respec- 
tive color maps. After extensive inspections, we resolve to 
chose a constant limit of mean slope (V606 — ^8i4)( r ) > 0.08 
mags /kpc 2 as a fiducial limit, that produces the cleanest sep- 
aration (Fig. [6] Using this limit we proceed and compute the 
fraction of active to total spheroids as a function of redshift 
for all the ellipticals in the Tadpole field up to z ~ 1.2. In or- 
der to avoid possible uncertainties arising from normal Pois- 
son variations and the use of photometric redshifts, we opt 
for binning out the data using redshift bins of size Az = 0.25, 
wide enough to smooth out major noise variations. In figure0 
we show the fraction of active to total systems as a function 
of their photometric redshifts. In order to incorporate the er- 
ror in the measurement of the photometric redshifts, 1 a error 
bars in figure were estimated by carrying out Monte Carlo 



simulations over the distributions. To estimate the uncertain- 
ties we recompute the histogram A times (N ~ 1000) in which 
the redshift of observation of the galaxies is modified such as 
z' u = Zo ± o:5z, in which a is a randomly generated number 
between [— 1,+1] and 8z is the nominal error estimation of 
BPZ equal to 8z = 0.15(1 +z). 

From figure0we report that the fraction of active spheroids 
is compatible with constant as a function of redshift (within 
the error bars). It is important to notice that the estimate for 
the lowest redshift bin, is quite uncertain as indicated by the 
large error bars, as th e Tadpole sample is incomplete for lower 
redshift galaxies (see lBeiutez et alJl2004l) . Regarding the ac- 
tual amount of the fraction of active spheroidals, we recognize 
that its absolute value is rather dependent in the limit used for 
selection. However, for the current limit adopted, this number 
mostly fluctuates around 30% — 40% in goo d agreement with 
the values obtained for the H DF-N for lKodama et alJ (119991) 
and Menanteau et al. (2001a). It is also worth noting that from 
tests performed on the sample, the shape of fraction distri- 
bution remains largely unchanged when adopting a different 
limit. 

It is appealing to analyze our estimate of the number of 
active spheroids in the context of their evolutionary history 
using the simple methodology introduced in this paper. The 
relatively constant to small rise in the fraction of active sys- 
tems is suggestive of the presence of a population of active 
spheroids, continuous at least up to z ~ 1 . This can also lead to 
the interpretation that spheroids are being constantly assem- 
bled since z ~ 1 . At the same time, the presence of spheroids 
with homogeneous colors and smooth color gradients is clear 
evidence of the presence of a relatively dominant population 
of normal evolved spheroids at all redshifts. 

It is possible to ask then, what type of objects are the blue 
spheroids that we detect? And, what would be their low red- 



Resolving spheroid galaxies from HST/ACS WFC UGC10214 ERO observations 



9 



FIG. 7. — The fraction of active systems to total as a function of redshift for 
spheroids in the Tadpole field. Error bars represent 1 o~ estimates obtained us- 
ing Monte Carlo simulations. The value for the first redshift bin is inprecise, 
as the sample is incomplete for low redshift galaxies. 



shift counterparts? A possible explanation, might be that we 
are seeing these galaxies as they enter into the E/SO class, 
possibly within the first Gyr, since their last star-formation 
activity. From the modeling of stellar synthesis populations 
jMenanteau et alJl2001albl) . it is believed that the existence 
of such populations is not long-lasting. We might be see- 
ing the predecessors or proto-ellipticals, in the last period of 
star-formation before becoming normal ellipticals as we know 
them at low redshifts and in clusters of galaxies. 

6. CONCLUSIONS 

In this paper, we have used newly available data from 
the Advanced Camera for Surveys to study the internal 
color properties of sample of 116 morphologically selected 
spheroids. This represents the first deep dataset of spheroids 
available from the ACS, and we used it to investigate its su- 
perb ability to resolve distant galaxies. 



Using an independent sample of ACS spheroids we have 
confirmed the presence of a population of spheroids with 
color inhomogeneities similar to those found in the HDFs. 
Using a model independent approach, we have introduced a 
new statistic, the mean slope (V(,o6 — hu)(f) which probes 
successfully the presence of blue cores and internal color 
inhomogeneities in spheroids. We compare our measure- 
ments with the HDF North with consistent results. Assum- 
ing that strong variations in the internal color dispersion of 
spheroids are linked to recent episodes of star formation, we 
use the mean slope to separate active and passively evolved 
systems as a function of redshift, and based on this, we esti- 
mate the fraction of active systems versus redshift. We found 
that within the uncertainties of our measurements, the frac- 
tion of active galaxies can be described as constant, with a 
trend to grow with redshift, with nominal values between 
~ 30% — 40% of the sample, consistent with previous re- 
sults. We take this as evidence for the continuous formation 
of spheroids since z ~ 1, while the data also shows the pres- 
ence of a population of old passively evolved ellipticals at all 
redshifts. 

An important part of the challenge that lies ahead is try to 
understand the significance of the blue spheroids and whether 
or not these represent the early stages of what we expect to 
be an elliptical as defined at z ~ 0. Additionally, alternative 
theories of galaxy formation need to be explored which can 
explain at the same time the presence of passive evolving sys- 
tems and recent star-forming spheroids. It also remains to 
be learned what are the detailed internal properties of ellipti- 
cals in clusters are from upcoming ACS observations at z ~ 1 
(Postman et al. 2004). 
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FIG. A8. — The distribution of synthetic PSFs created using Tiny Tim as function of the position on the geometrically corrected mosaiced ACSAVFC field. 
The plus symbols show the position of the synthetic PSFs in the WFC field. 



Menanteau, F., Ellis, R. S., Abraham, R. G., Barger, A. J., & Cowie, L. L. 

1999, MNRAS, 309, 208 
Menanteau, F, Jimenez, R., & Matteucci, F. 2001b, ApJ, 562, L23 
Papovich, C, Giavalisco, M., Dickinson, M., Conselice, C. J., & Ferguson, 

H. C. 2003, ApJ, 598, 827 
Peebles, P. J. E. 2002, in ASP Conf. Ser. 283: A New Era in Cosmology, ed. 

N. Metcalfe & T. Shanks 
Petrosian, V. 1976, ApJ, 209, LI 
Sandage, A. & Visvanathan, N. 1978, ApJ, 225, 742 

Sirianni, M., Ford, H. C, Illingworth, G. D., Clampin, M., Hartig, G., Becker, 
R. H., White, R. L., Bartko, F, Bemtez, N., Blakeslee, J. P., Bouwens, R., 
Broadhurst, T. J., Brown, R., Burrows, C, Cheng, E., Cross, N., Feldman, 
P. D., Franx, M., Golimowski, D. A., Gronwall, C, Infante, L., Kimble, 
R. A., Krist, J., Lesser, M., Magee, D., Mattel, A. R., McCann, W. J., 
Meurer, G. R., Miley, G., Postman, M., Rosati, P., Sparks, W. B., & 
Tsvetanov, Z. 2004, in preparation 

Tamura, N., Kobayashi, C, Arimoto, N, Kodama, T, & Ohta, K. 2000, AJ, 
119, 2134 



Trager, S. C, Faber, S. M., Worthey, G., & Gonzalez, J. J. 2000, AJ, 120, 165 
Tran, H. D., Sirianni, M., Ford, H. C, Illingworth, G. D., Clampin, M., 
Hartig, G., Becker, R. H., White, R. L., Bartko, F, Bemtez, N., Blakeslee, 
J. P., Bouwens, R., Broadhurst, T. J., Brown, R., Burrows, C, Cheng, E., 
Cross, N, Feldman, P. D., Franx, M., Golimowski, D. A., Gronwall, C, 
Infante, L., Kimble, R. A., Krist, J., Lesser, M., Magee, D., Martel, A. R., 
McCann, W. J., Meurer, G. R., Miley, G., Postman, M., Rosati, P., Sparks, 
W. B., & Tsvetanov, Z. 2003, ApJ, 585, 750 
Treu, T. & Stiavelli, M. 1999, ApJ, 524, L27 

Treu, T., Stiavelli, M., Casertano, S., M0ller, P., & Bertin, G. 1999, MNRAS, 
308, 1037 

van Dokkum, P. G. & Franx, M. 1996, MNRAS, 281, 985 
White, S. D. M. & Frenk, C. S. 1991, ApJ, 379, 52 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 
Zepf, S. E. 1997, Nature, 390, 377 



Resolving spheroid galaxies from HST/ACS WFC UGC10214 ERO observations 



11 



> 



CO 



CD - 
O d 



01 



10 



02 



03 



12 



04 



0.2 0.4 0.6 0.1 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.1 

Central Concentration 



0.2 0.4 0.6 0.1 



FIG. A9. — The internal color dispersion 8(V — I), used to measure color inhomogeneities, as a function of the central concentration, C. Open circles represent 
simulated spheroids to estimate the bias le vel i ntroduced by the PSF. Open squares observed ones. Each panel shows the results of artificial spheroids for each 
PSF positions using the description of Fig. lA8l 
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FIG. A10. — The mean slope (V606 — hu){>")i use d to trace central blue light, shown as a function of the central concentration, C. Open circles represent 
simulated spheroids to estimate the bias le vel i ntroduced by the PSF. Open squares observed ones. Each panel shows the results of artificial spheroids for each 
PSF positions using the description of Fig. lA8l 



12 F. Menanteau & The ACS Science Team 



Table B3. The catalog of selected spheroids from UGC102104 



ID 


RA(J2000) 


DEC(J2()00) 


F814W 


S(y-i) 


mean slope (V 6 06 — /gi4)('") 


BPZ 


A 


C 


1041 


1 6-06-16 49 


+SS-23-44 19 


21.877 


n 0691 +0 0019 


n 0587 ±0 001? 

U.V/JO / 1 \7 . \J\J ± _ 


0.23 


0.054 


0.538 


1 076 
1\J /V 


i ft-oft oi ss 


-+.SS-94-S4 on 

TJJ .ZH-. J4.UU 


22 982 


01 4- 0079 
VJ.Vj ItO In U.UU / a 


09Q9 4- 0^"*S 


67 




392 


1 1 oi 


1 ft*Oft*OQ Rft 


-l.SS-94-98 ftn 


21 767 


0^0^ -1- 009ft 


04RS 4- 009ft 


21 




o iss 

U.JOJ 


1 109 


1 U . UU .17.10 


-4-SS-94.-17 ft^ 

TJJ.il.l / .UJ 


23 713 


1001 4-0 09RS 
\j. io"i zn u.uioj 


o 1 si s 4- o o^so 

\j . i j i j m u.ujju 


47 


067 


275 


1142 


16:06:06.01 


+55:24:53.35 


23.408 


0.1284 ±0.0277 


-0.0982 ±0.0386 


0.41 


0.047 


0.365 


1157 


16:06:13.95 


+55:24:57.49 


21.760 


0.0964 ±0.0086 


-0.0765 ±0.0030 


0.78 


0.052 


0.466 


1200 


16:06:04.61 


+55:25:58.07 


23.408 


0.5834 ±0.0079 


-0.2479 ±0.0304 


0.85 


0.051 


0.477 


1207 


16:06:01.80 


+55:26:17.99 


20.447 


0.1038 ±0.0013 


-0.0751 ±0.0009 


0.71 


0.081 


0.653 


1266 


16:06:18.02 


+55:24:31.97 


23.362 


0.0802 ±0.0219 


0.1880±0.0419 


0.61 


0.004 


0.451 


1275 


16:06:19.02 


+55:24:25.81 


22.822 


0.0882 ±0.0080 


-0.0050 ±0.0077 


0.74 


0.041 


0.425 


1282 


16:06:15.43 


+55:24:49.85 


23.337 


0.0563±0.0151 


0.1083±0.0159 


0.83 


0.130 


0.476 


1295 


16:06:03.27 


+55:24:57.76 


20.654 


0.0701 ±0.0058 


-0.0658 ±0.0473 


0.78 


0.090 


0.428 


1364 


16:06:15.43 


+55:24:52.52 


23.303 


0.0330 ±0.0469 


0.1256 ±0.0582 


0.78 


0.128 


0.702 


1475 


16:06:19.52 


+55:24:30.48 


23.856 


0.0490 ±0.0149 


0.1673±0.()175 


0.67 


0.040 


0.379 



Complete table available at http: //acs .pha. jhu. edu/~f elipe/e-prints 



APPENDIX 

A. THE EFFECT OF THE ACS/WFC POINT-SPREAD FUNCTION 

It is known that systematic variations present in the HST PSF may influence studies which depend on small aperture photometry 
or small galaxies with bright cores. Our main concern is the PSF variation as a function of wavelength, which might induce 
spurious centrally concentrated color inhomogeneities on sharply peaked light profiles. To investigate the bias level introduced 
by variations of the PSF in our i nternal color analysis of sp heroids, we carried out Monte Carlo simulations following a procedure 
similar to the one described in Menanteau et al. (2001a). This consists of creating artificial galaxies using the IRAF package 
ARTDATA resembling the sizes and brightness as of those in our study. These were subsequently analyzed in the same fashion 
as the observed data. Initially we attempted to create artificial de Vaucouleurs profiles using a combination of averaged observed 
bright stars in the background of the Tadpole field as PSF. However, it was soon evident that there was a systematic variation in 
the PSF FWHM as a function of position in the WFC field. In addition, the small number of bright stars in The Tadpole field 
made the use of stars as real PSFs too challenging. The PSF variation across the WFC field has been extensively documented by 
iRristl (120031) . Although the PSF variat ions in ACS are less than other HST cameras, its does varies over the field, most notably 
for the WFC (see Fig. 1 in Krist 2003). Therefore, we chose to utilize the synthetic PSFs generated with the Tiny Tim software 
as it does take into account field variations. 

We created a grid of 4 x 4 PSFs at evenly distributed positions across both WFC chips with Tiny Tim for the V606 and I%u 
bands. These were incorporated into ACS mock raw images and later processed through Apsis using the same parameters as 
for the observed data. These yield a geometrically corrected grid image of synthetic PSFs. In Fig. lA8l we show a diagram with 
the final positions. Using this grid, spheroids were synthesized at each PSF location with de Vaucouleurs profiles using IRAF's 
ARTDATA. A set of simulated spheroids were generated for the Vgoe and 7gi4 bands trying to mimic as much as possible the 
observed properties of spheroids. These include a range of half-light radius (r e ) from 0.1" — 0.7" which correspond to a physical 
length of - 0.75 - 5.5 kpc (H = 71 km s^Mpc" 1 , £2 = 0.27, Q. A = 0.73) and a range of magnitudes consisted with our 
sample (18 < 7x14 < 24). We processed the synthetic spheroids using the exact same procedure employed for real spheroids and 
calculated both estimators, 8(V — I) and the mean slope (Ve06 — ^814) ( r ) f° r the whole sample. To simplify the simulation we 
compute the mean slope as function of the object apparent size. As the influence of the PSF in the resolved colors is expected to 
be a strong function of the 'peakiness' of light distribution o f the galaxy, we probe the influence of the PSF on real and simulated 
galaxies against the central concentration parameter (sec tion|23). 

The results of this exercise are shown in Figs. IA9I and lAlOl where we plot the values of 8(V —I) and the mean slope 
(V606—hu)(r) as a function of the central concentration respectively for the 16 PSF pointings. In all 16 positions we also 
plot the observ ed da ta (open squared) to compare with the simulation results. When confronting real and simulated values of 
8(V — I) in Fig. lA9l we can see that over the observed range of C for the ACS spheroids (i.e. 0.3 < C < 0.7) the values of 8(V — I) 
for the simulated galaxies as significantly lower that those obtained for the observed spheroids. Only for values of C > 0.8 we 
enter into the regime in which the PSF can affect the recovered values of 8 (V — I) as the PSF becomes important. However, our 
spheroids do not lie in this range. We notice some variation depending on the PSF position used, but this effect is too small for 
the scale of the real valu es of 8(V —I). We also compare the central concentration with the mean slope (V606 — ^814) (r) using the 
same symbols as of Fig. lA"9l The simulated galaxies have mostly a uniform distribution of values of mean slope (V606 — hu)( r ) 
as a function of C, also lower than the observed for the real data. We conclude that the bias level introduce by the variations in 
the WFC PSF will not seriously affect any of our conclusions. 
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